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SOME EFFECTS OF NONLINEAR VARIATION IN THE DIRECTIONAL-STABILITY 
AND DAMPING-IN- YAWING DERIVATIVES ON THE 
LATERAL STABILITY OF AN AIRPLANE 1 


By Leoxasd Stsbxfield 


SUMMARY 

A theoretical investigation has been made to determine the 
effect of nonlinear stability derivatives on the lateral stability 
of an airplane. Motions were calculated on the assumption 
that the directional-stability and the damping-in-yawing 
derivatives are functions of the angle of sideslip. The. applica- 
tion of the Laplace transform to the calculation of an airplane 
motion when certain types of nonlinear derivatives are present is 
described in detail. The types of nonlinearities assumed 
correspond to the condition in wh ich the values of the directional- 
stability and damping-in-yawing derivatives are zero for small 
angles of sideslip. 

The results of the investigation indicated that under certain 
conditions the nonlinear stability derivatives assumed in the 
analysis caused a motion which had different rates of damping 
for the large and small amplitudes of motion, with very little 
damping at the small amplitudes. In general, the period of the 
resultant oscillation increased with time. 

INTRODUCTION 

Recent flight tests of several airplanes designed for high- 
speed high-altitude flight have indicated neutrally damped 
lateral oscillations of small amplitude generally referred to 
as snaking. Upon examination of the flight records, the 
decrement of the oscillatory motion is found in some cases 
to be different for the large and small amplitudes of motion 
with, a neutrally stable oscillation occurring at the small 
amplitudes. One of the explanations offered for the cause 
of this type of motion is that some of the stability derivatives 
are nonlinear: that is, the derivatives have different values 
for the large and small amplitudes of motion. The non- 
linearity could be caused by boundary-layer effects or flow 
separation due to poor fairing at the junction of the. tail 
surfaces. 

The present report represents a preliminary investigation 
of the effect of the presence of two nonlinear stability 
derivatives, the directional-stability derivative 0^ and the 
damping-in-yawing derivative C Ur , on the motion of an air- 
plane. These derivatives were selected for the analysis 
since the damping of the oscillation is a function of C %r and 


since C* r depends upon the C„p contributed by the tail. 
The derivatives G„. and C n , were both assumed to be func- 
tions of the sideslip angle /3. Calculations were made of the 
airplane motion due to a disturbance in sideslip for three 
different types of nonlinearities. 
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angle of roll, radians 
angle of yaw, radians 

angle of sideslip, radians except where noted 
in figures ( vjV) 

yawing angular velocity, radians per second 

(.dt/dt) 

rolling angular velocity, radians per second 
(d<t>!dt) 

sideslip velocity along the F-axis, feet per 
second 

airspeed, feet per second 

mass density of air, slugs per cubic foot 

dynamic pressure, pounds per square foot 

S'**) 

wing span, feet 
wing area, square feet 
weight of airplane, pounds 
mass of airplane, slugs (Wjg) 
acceleration due to gravity, feet per second per 
second 

relative density factor (mjpSb) 
inclination of principal longitudinal axis of 
airplane with respect to flight path, positive 
when principal axis is above flight path at 
the nose, degrees 

angle of flight path to horizontal axis, positive 
in climb, degrees 

radius of gyration in roll about principal 
longitudinal axis, feet 
radius of gyration in yaw about principal 
vertical axis, feet 

nondimensional radius of gyration in roll about 
principal longitudinal axis ( fxfb ) 


1 Supersedes XACA TX 2233, "Some Effects of Xoniinear Vadatton hi the Directional-Stability and Damping-in-Yawing Derivatives on the Lateral Stability of an Airplane” by 
Leonard Stemfleld, 1060. 
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K z , 
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nondimensional radius of gyration in yaw 
about principal vertical axis ( k z Jb ) 
nondimensional radius of gyration in roll 
about longitudinal stability axis 
(V-Ky, 2 cos 2 v+Kz* sin 2 v) 
nondimensional radius of gyration in yaw 
about vertical stability axis 
(V#V cos^+Kx/sIES) 
nondimensional product-of-inertia parameter 
((Kz^—KxJ) sin tj cos ij) 


trim lift coefficient ~ ~ ^ 

. - . , /Rolling moment\ 

rollmg-moment coefficient ( J 

. /Yawing moment\ 
yawing-moment coefficient ^ qSb ) 

lateral-force coefficient ( LateraHorc e.^ 



C n 


t 


*0 


n> . 


yawing-moment constant 
time, seconds . 

nondimensional time parameter based on span 
(Yt/b) 


differential operator 



o- operator in Laplace, transformation 

T h time for amplitude of oscillation to damp to 

one-half its original value, seconds 

The subscript 0 is used to indicate initial conditions and 
u bar is used to denote variables in the operational equations. 


ANALYSIS 

NONLINEAR STABILITY DERIVATIVES 

The assumptions made with regard to the nonlinearity of 
the stability derivatives C„» and C„ r are shown in figure 1. 
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(a) Case 1. (b) Case 2. (e) Case A 

FinuBE 1 — Three types of nonlinear stability dc-lvatlves assumed In Urn analysis. 


For all three cases presented in the figure, C\ is equal to zero 
for — 2°</S<C2°, a region which is subsequently referred to 
as a dead spot. Thus when the airplane is within the dead 
spot, the value of the directional stability derivative C Uf 
is zero. Since the damping-in-yawing derivative (\ r is a 
direct function of C nf contributed by the tail, C, r was also 
assumed to be zero for values of — 2°<j8<C2°. In the region 
outside of the dead spot, each one of the cases represents a 
different type of variation of C n with /3 in order to simulate 
the effect of several possible flow conditions on the side force 
acting on the vertical surface. For cases 1 and 2, C„ t = 0.28 
and for case 3, C»,=0.41. The corresponding value of 
Cn T for all three cases is —0.39. It should be noted in figure 1 
that for cases 2 and 3, C»=0 at of 2° and —2°, whereas 
for case 1, C„ has a finite value at 0 of 2° and —2°. 


METHOD OF CALCULATING MOTION 

Since the nonlinearities shown in figure 1 can be treated 
as linear derivatives of different values within and outside of 
the dead spot, the airplane motion is calculated on the basis 
of classical linear theory. The equations of motion and the 
general method of calculating the motion of an airplane are 
given in references 1 and 2. The methods of references 1 
and 2 are based on the Laplace transformation which in- 
herently takes into account the initial conditions of the 
problem. Because the Laplace transformation considers 
the initial displacements and initial velocities of the problem, 
this method is directly applicable to the calculation of the 
motion of an airplane which has nonlinear derivatives 
similar to the derivatives presented in figure 1. 

The nondimensional linearized lateral equations of mo- 
tion, referred to the stability axes, are for rolling, for yaw- 
ing, and for sideslipping, respectively: 

2^K x Wf<p+K xz DM ) = <?,/ +-* C tf D>4>+ \ (\ 
2,^K z WM+Kx Z DM)= 1M+\CnDrf-Y C. e ^ 

D t> -f)=C Y ^-\-- I ^Cy v Di > 4>-\- C L <t>-\- Y' r , 

(C t tan y)t 

(0 




SOME EFFECTS OF NONLINEAR VARIATION IN C Kjt ANT) C„ r ON LATERAL STABILITY OF AN AIRPLANE 


1011 


The Laplace transformation of equations (1) , with, tlie use of the symbol <r for the operator, is 
( S/IalvyV ' — G — gC'l r <^\p — Cj^jS = 2jJ.l ) K X 2 [ff $0+(LW)o] — g- Cf p <£ 0 + 2 flfiKxZ W'Po J r(Db 1 f)o] — ■ 

^ C» p cJ<p + ( 2 ^jAz 2 tr s — — ( 7 Sr (r^ — O x J5=2ni ) K xz [e , $o - f - (i^ 6 $)oI — ^Cx II 4>a~\-2(tb&-z 2 { a ’4'a J r(Dt,ip) 0 \ — ypn r } ! r v \ — “ 
( — ’TyCr^o’ — CL^<p-\-(2fi 6 <r — ^-Cr r cr — Cl tan T^+(2^scr— Cr jj )jS= — 9 ^V p <£o+^ 2^ — ^Cr r ^)&>+ 2^400 


K 2 > 


Equations (2) represent three simultaneous algebraic equations which can be solved for ft, <j>, f, and their derivatives by 
the method of determinants. For example, 


■* Al /(<r ) 

where A is the characteristic lateral-stability equation 

(^rltr' 4 — L tr 3 — f— Ccr 5 -\-D<t-\-E)g 

Ai = j.'iiO" 1- l _ -Bicr 5 -(-Ci(r i -\~Dy<r-\-Ei 


( 3 ) 


and 


The expressions for H, B, C, D, and E, in terms of the mass and aerodynamic parameters of the airplane, are given on 
pages 27 and 28 of reference 1. The coefficients of the A, equation are 

Ai=p 0 [S^iK^KZ—Kx^] 

Bi = <t>a [4 rfCdKSKzi-Kxz*)] +{TM)« [2rfOr p (K x i K z t -K xz *)] + 4>a [4m» 2 C £ tau y (K X 2 K Z *—K XZ *)] + 

(ZV)o [2^\Cy r -4^) (Kx 2 K/—K xz i )} + p 0 [2 rf(K xz C lr -K x *0 Xr ) +2rf(K xz C Xv -KfO t j\ 

c I =0 o [M6C i (ic xz C: r -ii:x 2 c„ r )+M S c t (£:x Z c % --^/c f j]+(D^) o [|^x s (c % c rr -c tt c i - p )+ 

2 l x b i K x \2K^C L -C Xs ) +2rfK xz (C lp -2K xz C L ) +\n>K xz (C t C Yp -C,Cr r )J+h [y- b C L tan y(K xz C.\-O x K x *) + . 

MiCitail y (K XZ C a p — FL Z C j p )] +(-Dj^)o Pi>K xz (C Xp C y r — C n C r p ) +;y t r Cy p — C Ip C r r ) + 

4/t» 2 C t tan y (K x ~K z i — Kx^) -f- 2 \L b (E. z C [p — K XZ C „ p ) Po Jjy (Cj p C„ r — C iC „ p ) J-f- 

Cn t [ — 1*I>K XZ Cy p ~\~ V-bEix'C Y r — 4fij-i£x ? [ 

Ei = < f>o jjj C l (C i p C — Cj r C Kp ) J+fCWMp&CL tan y(K x 2 C Xp — K xz Ci p ) + v-tP l(K xz C i r — K X ~C + 

'Pn [| C h t^ny (C t C n -C l C Xp )^{B^U^C^{E z -C l -K xz C Xr )+^C L t&iiy{K xz C np -KFC lp )} + 

C „ c |^— Cr p Cf r — 2fi b K xz C fi b E x i C l tan y +MtCi p — ^C t p C j- r J 

E,=C Wc C L C lr -± C lp C L tan 7 ) 

The solution of equation (3), which will result in a time history of /3 as a function of s t . is obtained from the Heaviside 
expansion theorem (reference 3) : 

g^v /W x,,, m 

where A* are the roots of F(a-) set equal to zero. Similar solutions are derived for <j>, f, D d <t>, Dt,p, and D b j 3. The time 

scale is readily converted from s b units to t units by the equation t=^iv 

The values of the stability derivatives and mass characteristics used in the calculations are presented in table I. The 
table is divided into two columns which differ only in the values of C Xf and C„ r of the airplane for the cases where the 
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TABLE I 

STABILITY DERIVATIVES AND MASS CHARACTERISTICS 
OF THE AIRPLANE CONSIDERED IX THE ANALYSIS 


Derivative or characteristic 

Outside of 
dead spot 

Within dead 
spot 

T17S, 

80 

SO 

"t - — 

101.1 
0. 01X189 

10L1 

0.00089 

v, ft/a>c . 

753 

753 

Ci-.'..'. . 

h t ft 

0.318 

27.7 

0.318 

27.7 

r.dfS 

0 

0.0573 

0 

0.0573 

jKV ... .. 

0.0009 

a 0060 


— 0. 462 

-0 . 402 


-a 0155 

-0.0155 


—0. 126 

-a 126 


0 

0 

Cr , per radian 

v, 

, pur r itfl lari 

• 0 

-2.-0, 0,2.0 
-a 392 
0.28 

0 

-2. 0, 0, 2. 0 
0 
0 

Cup (case 3), pcrradlan . 

an 

0 




airplane is either outside of or within the dead spot. From 
the analytical solution of the motion, based on the mass and 
aerodynamic characteristics of the first column of table X 
and an initial condition of /S=5°, the time history of 0 was 
computed for several values of s b until the value of s b for 
which /3=2° was reached. For values of greater than the 
s b which results in /3=2°, this analytical solution is incorrect 
since the airplane has now entered into the dead spot and the 
values of C nft and are zero. Thus, a new solution must 
be calculated with the use of the values given in the second 
column of table I with new initial conditions. The new 
initial conditions are determined by substituting the value of 

at which d=2° in the original analytical solutions of <fi, 
i P, Db4>, Drf, and D fj p. Once these initial conditions are 
known, another set of analytical solutions are computed for 
/3, <f>, and their derivatives from equations (3) and (4). 
This procedure is followed every time /3 crosses through 2° 
or —2°. The final resultant motion in sideslip is the sum of 
all the analytical solutions in /3, each one of which is correct 
only for a particular interval of time. 

The constant C„ e is introduced into the yawing-moment 
equation of equations (1), since the value of the yawing- 
moment coefficient due to sideslip is +C’„ c for the con- 
dition of the airplane having the dead spot in cases 2 and 
3 of figure 1. The values of C B( . are [0.00977 { and |0.0143j 
for cases 2 and 3, respectively. The sign of C„ c is opposite 
to that of (3. For case 1 of figure 1, C* c =0. 

It is apparent that the procedure employed is a time- 
consuming process and subject to the possibility of many 
computational errors due to the magnitude of the computa- 
tions. The final solution can be obtained, however, in a 
relatively short time through the use of automatic- digital 
computing machines. 


RESULTS AND DISCUSSION 

The effect of the nonlinear stability derivatives on the 
lateral motion was investigated for the airplane described 
by the" mass an<f aerodynamic characteristics given in table 
I, with three different values for the damping of the lateral 
oscillation as calculated on the basis of derivatives constant 
with amplitude. Since the damping was varied arbitrarily 
by assuming different values for the angle of inclination of 
the principal longitudinal axis of the airplane to the flight 
path r), three values of i?, —2°, 0°, and 2°, were selected 
which correspond to a damping of the lateral oscillation, 
expressed in terms of T m , of 5.6, 3.0, and 1.8, respectively. 
The motion of the airplane in sideslip, due to an initial 
disturbance in sideslip of 5°, for the three values of y is 
shown in figure 2. Since these motions are calculated on 
the assumption of derivatives constant with amplitude, the 
amplitudes of the motion decrease exponentially with time 
and wilL eventually reduce to zero. As can be noted in the 
first column of table I, the C*. for cases 1 and 2 is 0.28; 
whereas the C». for case 3 is 0.41. The motions presented 
in figure 2 are for C*.=0.28; however, the motions for 
C^=0.41 would exhibit oscillations of approximately the same 
damping and a slightly smaller period. 

The motions of the airplane in sideslip, showing the effect 
of the nonlinearities illustrated in figures 1(a), 1(b), and 1(c), 
are presented in figures 3 to 5, respectively. In all cases, 
an initial disturbance in sideslip of 5° was assumed. The 
pronounced effect of the nonlinearities on the lateral motion 
is noted by a comparison of figure 2 and either one of figures 
3, 4, or 5. In all three figures (figs. 3 to 5) the motion for 
v=2°, the most stable case, approaches a constunt value. 
The analytical solution of the motion for the case of rj=2° 
in figure 3 indicates that, within the dead spot, the airplane 
will oscillate at a period of 6.56 seconds and 7’i / i=3.38 
seconds and will eventually approach the value of 0= — 0.0092°. 
Similar motions would be obtained for the ease of rj=*2° in 
figures 4 and 5. As 17 is decreased, the damping of the 
oscillatory motion depends upon the nonlinearity assumed 
and the values of n. In figure 3, the motion for 17=0° 
damps at a slow rate at the large amplitudes until the oscil- 
lation reaches an amplitude of approximately 2.4° where 
the damping of the oscillation is zero. The period of the 
oscillation increases from L5 to 1 .85 seconds. For the case of 
t)=— 2^-a very lightly damped oscillation is apparent within 
the first few seconds and the airplane may be considered to 
be neutrally stable at an amplitude of F4.5 0 . In figures 
4 and 5 the motion for 17=0° clearly indicates that the damp- 
ing is decreasing as the amplitude decreases and the period 
of the oscillation increases; for t )= — 2°, the oscillatory motion 
is slightly unstable. A neutrally stable oscillation would be 
expected to occur in figures 4 and 5 for the combinations of 
a value of rj between 0° and —2° and the dead spot assumed 
in the calculations or for 17=— 2° and a smaller dead spot. 
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In general, the results indicate that the damping of the 
lateral oscillation calculated with the use of derivatives con- 
stant with amplitude is a determining factor in the type of 
motion obtained where nonlinear derivatives are present. 
As the inherent damping of the lateral oscillation decreases, 
a smaller dead spot will result in a neutrally stable oscillation. 
Obviously, if the inherent damping is zero, a neutrally stable 
oscillation already exists with zero dead spot. 


Some. additional calculations w F ere made for the case where 
the airplane is disturbed within the dead spot. The motions 
for an initial condition of /3 = 1° were computed for ij =— 2° 
and 0° with the assumption of the nonlinearity described in 
figure lib). The results are presented in figure 6. It should 
be noted that the only difference between figures 4 and 6 is 
the initial condition assumed in the calculations. In figure 6, 
the motion for 77 = — 2° is unstable and gradually approaches 
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Fracas 6.— The eflect of the nonlinear derivatives described In figure 1(b) on the motion of an airplane. Initial dtstnrbance In sideslip of 1 ". 


the amplitude and period of the motion for the case of 17 = — 2 ° 
in figure 4. The motion for 17 = 0 ° in figure 6 is slightly un- 
stable and will probably increase until its amplitude and 
period are in close agreement with the motion for the case 
of ij= 0 ° in figure 4. Calculations have indicated that the 
oscillatory motion of the airplane within the dead spot will 
double amplitude about every 4 seconds for 17 = — 2° and 
about every 30 seconds for 17 = 0 °. If the motion is unstable 
within the dead spot, either the airplane motion will be neu- 
trally stable with an amplitude equal to or greater than the 
amplitude of the dead spot or the motion will be unstable. 

The loss in damping and the increase in period which ap- 
peared in some of the lateral oscillations in figures 3 to o can 
be attributed to the type of nonlinearity assumed. From 
classical dynamic stability theory, it is well known that the 
damping of the oscillation is a function of and the period 
of the oscillation is a function of C np . If the airplane is con- 
sidered as a mass-spring dashpot system, O n$ is the equiv- 
alent spring constant of the system and C nr corresponds to 
the damping constant contributed by the dashpot. Thus as 
C n p is reduced the period increases, and as £7„ r is reduced the 
damping decreases. 


CONCLUDING REMARKS 

The results of the investigation made to % determine the 
effect of nonlinearities assumed in the analysis on the lateral 
stability indicate that under certain conditions a motion is 
obtained which has different rates of damping for the large 
and small amplitudes of motion, with very little damping at 
the small amplitudes. In general, the period of the resultant 
oscillation increases with time. 

/ 
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